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We study quantum interference effects on the transition strength for strongly localized electrons 
hopping on 2D square and 3D cubic lattices in a magnetic field B. In 2D, we obtain closed-form 
expressions for the tunneling probability between two arbitrary sites by exactly summing the cor- 
responding phase factors of all directed paths connecting them. An analytic expression for the 
magnetoconductance, as an explicit function of the magnetic flux, is derived. In the experimen- 
tally important 3D case, we show how the interference patterns and the small-B behavior of the 
magnetoconductance vary according to the orientation of B. 
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Quantum interference (QI) effects between different 
electron paths in disordered electron systems have been 
a subject of intense study Jl|-^] because they play an 
important role in quantum transport; for instance, the 
QI of closed paths is central to weafc-localization phe- 
nomena 0. Recently, a growing interest exists on the 
effects of a magnetic field on strongly localized elec- 
trons with variable-range hopping (VRH) where striking 
QI phenomena have been observed in mesoscopic and 
macroscopic insulating materials. This strongly local- 
ized regime J2j-|| is less well-understood than the weak- 
localization case. Deep in the insulating regime, the ma- 
jor mechanism for transport is thermally activated hop- 
ping between the localized sites. In the VRH regime, 
localized electrons hop a long distance (the lower the 
temperature is, the further away the electron tunnels) 
in order to find a localized site of close energy. The 
conductance of the sample is governed by one critical 
phonon-assisted hopping event 0. During this critical 
tunneling process, the electron traverses many other im- 
purities since the hopping length is typically many times 
larger than the localization length. It is important to 
emphasize that the electron preserves its phase memory 
while encountering these intermediate scatterers. This 
elastic multiple-scattering is the origin of the QI effects 
associated with a single hopping event between the initial 
and final sites. The tunneling probability of one distant 
hop is therefore determined by the interference of many 
electron paths between the initial and final sites ^ M . 

In this paper we investigate the QI of strongly localized 
electrons by doing exact summations over all directed 
paths between two arbitrary sites. For electrons prop- 
agating on a square lattice under a uniform potential, 
we derive an exact closed-form expression for the sum- 
over-paths. We also obtain an explicit formula for an 
experimentally important case, much less studied theo- 
retically so far: the interference between paths on a 3D 



cubic lattice. In the presence of impurities, by computing 
the moments of the tunneling probability and employing 
the replica method, we derive a compact analytic result 
for the magnetoconductance (MC), which is applicable in 
any dimension. Our explicit field-dependent expression 
for the MC provides a precise description of the MC, 
including the low and high field limits. The period of os- 
cillation of the MC is found to be equal to hc/2e. Also, 
a positive MC is clearly observed when turning on the 
field B. When the strength of B reaches a certain value, 
B sa , which is inversely proportional to twice the hop- 
ping length, the value of the MC becomes saturated. At 
very small fields, for two sites diagonally separated a dis- 
tance r, the MC behaves as: rB for quasi-lD systems, 
r 3 / 2 B in 2D with B = (0,0,5), and rB [r 3 / 2 B] in 3D 
with B parallel [perpendicular] to the (1, 1, 1) direction. 
The general expressions presented here (i) contain, as 
particular cases, several QI results derived during 

the past decade (often by using either numerical or ap- 
proximate methods), (ii) include QI to arbitrary points 
(m, n), instead of only diagonal sites (m, m), (iii) focus 
on 2D and 3D lattices, and (iv) can be extended to also 
include backward excursions (e.g., side windings) in the 
directed paths. 

Exact results in this class of directed-paths problems 
are valuable; and, for instance, can be useful to study 
other systems: (1) directed polymers in a disordered sub- 
strate, (2) interfaces in 2D, (3) light propagation in ran- 
dom media, and (4) charged bosons in ID. 

To study the magnetic field effects on the tunnel- 
ing probability of strongly localized electrons, we start 
from the tight-binding Hamiltonian H — W^2 i c\ci + 

VJ2(ij) c\ c i eiMi > where V / W < 1 and A v = 2lr Ji A ' dl 
is 27r times the line integral of the vector potential along 
the bond from i to j in units of $o = hc/e. Con- 
sider two states, localized at sites i and / which are 
r bonds apart, and the shortest-length paths (with no 
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backward excursions) connecting them; i.e., the directed- 
path model. By using a locator expansion, the Green's 
function (transition amplitude) between these two states 
can be expressed as T if = W{V/W) r S {r \ with = 
^ r exp(i$r), where T runs over all directed paths of 
r steps and $r is the sum over phases of the bonds 
on the path. This directed-path model provides an 
excellent approximation to Tif in the extremely local- 
ized regime |^-^] since higher-order contributions involve 
terms propotional to W(V/W) r+21 (I > 1), which are 
negligible because (V/W) 2 is very small. Quantum in- 
terference, contained in S^ r \ arises because the phase 
factors of different paths connecting the initial and final 
sites interfere with each other. In this work, we focus on: 
(i) the computation of the essential QI quantity S^ r \ (ii) 
the derivation of the MC in the disordered case, and (Hi) 
the study of the full behavior of the MC — including the 
scaling in the low-field limit and the occurrence of satu- 
ration. It is important to keep in mind that the effect of 
a magnetic field on the MC follows the behavior of S^- r \ 
Quantum Interference on a 2D square lattice. — Let us 
choose (0, 0) to be the initial site and focus on sites (to, n) 
with m, n > 0. For forward-scattering paths of r steps, 
which exclude backward excursions, ending sites (m, n) 
satisfy m + n — r. Let S m<n (—S^) be the sum over all 
directed paths of r steps on which an electron can hop 
from the origin to (to, n), each one weighted by its corre- 
sponding phase factor. Employing the symmetric gauge 
A = (—y,x)B/2, and denoting the flux through an el- 
ementary plaquette (with an area equal to the square 
of the average distance, which is typically equal to or 
larger than the localization length, between two impuri- 
ties) by <p/2w, it is straightforward to construct the recur- 
sion relation: S m , n = e~ m ^ 2 S m -i, n + e im< t>/ 2 S m>n - V 
The factors in front of the 5"s account for the pres- 
ence of the magnetic field. Enumerating the recur- 
sion relations for Sk ni n (k n = m — 1, ... ,0) successively 
and using So^ n = 1, we obtain the following relation 
Sm,n = EZ=o etk " H2 e-**™-*")"*/ 3 Sfc n ,„-i- This equa- 
tion states that the site (to, n) can be reached by moving 
one step upward from sites (k n , n — 1) with < k n < to, 
acquiring the phase ik n (j>/2; then traversing to — k n steps 
from (k n ,n) to (m,n), each step with a phase —in<p/2. 
By applying the above relation recursively and utiliz- 
ing S m fl = 1, S m ^ n for to, n > 1 can be written as 
S m , n ((p) = exp(-imn<j>/2) £ m ,„(0), and 

ww>)=E E e i ^ + - +fc - i+fc -)* (i) 

fc„=0fc„_i=0 fei=0 

Notice that each term in the summand corresponds to 
the overall phase factor associated with a directed path. 
When = 0, S m>n (0) = C r m = N is just the total 
number of r-stcp paths between (0,0) and (to, n). 

After some calculations we obtain one of our main re- 
sults, a very compact and elegant closed-form expression 
for S m . n ((f>): 



Notice that the symmetry S m ^ n — S n>m [apparent in 
Eq. (2)] is due to the square lattice geometry. In the 
very-low-flux limit <f> <C 1, the logarithm of S m ,n, calcu- 
lated exactly to order (ft 2 (and omitting In N) , is 

In S m ,n{4>) — ~^ mn ( m + n + 1) ( 3 ) 

and thus we obtain the familiar Q harmonic shrinkage 
of the wave function. 

Sm,m has the richest interference effects because the 
number of paths ending at (m, to) and the areas they 
enclose are both the largest. We therefore examine 
more closely the behavior of the quantities I^m (0) = 
Sm,m(<l>) = nr=i(sin^0)/(sin|0). I 2m {<j>) obeys 
the following properties: (i) 2it (47r) periodicity in for 
even (odd) m, (ii) 7 2 m(27r — <f>) = I 2m (0) for < cf) < ir 
with to even, and (iii) /2m(27r ± <p) — —l2m(4 > ) f° r 
< 4> < 2n with m odd. Furthermore, the zeros of 
J 2m (0) are given by <f> = 2irs/t, for ^ < i < ^, 
with < n < and the s's are prime to each allowed 

t. From the physical viewpoint, these flux values pro- 
duce the complete cancelation of all phase factors (i.e., 
fully destructive interference). Indeed, as the magnetic 
field is turned on, l2m(4>) rapidly drops to its first zero 
at <J)/2tt = 1/2to and then shows many small-magnitude 
fluctuations around zero. 

Effects of disorder. — To incorporate the effects of ran- 
dom impurities, we now replace the on-site energy part 
(first term in H) by X)i e « c I c 2- Now the are inde- 
pendent random variables which can take two values: 
+W with probability fi and —W with probability 
where fi + v = 1. Due to disorder, the transition am- 
plitude becomes = W(V/W) r J m<n , with J m ,n = 
Er[IIjer( — W/tj) ]e l * r . For all directed paths ending 
at (to, n), electrons traverse r sites (the initial site is ex- 
cluded). Each site visited now contributes an additional 
multiplicative factor of either +1 or —1 to the phase fac- 
tor. Therefore, for a given path T, the probability for 
obtaining ±e I$r is [(n + v) r ± (fi — v) r \/2. By exploiting 
Eq. (1), we derive a general expression for the disorder av- 
erage of the tunneling probability (i.e., the transmission 
rate) as ( | J m ,„(0)| 2 > - N+^-uf^S^cfy-N), where 
(■ ■ •) denotes averaging over all possible configurations of 
impurities. It is important to stress that the conductivity 
between i and / is proportional to ( | J m ,n{4>)\ 2 ) @ ||- 

For the most studied case so far: fi = v = 1/2, 
we can obtain analytical expressions for the moments 
( \J m ,n{ ( l ) )\ 2p ) for any value of p. In general, ( | J m ,„(</>)| 2p ) 
consist of terms involving N k with k = 1, . . . ,p. Here- 
after, we omit the subscripts in J and S, and focus on 
the leading terms (oc N p ), since they provide the most 
significant contribution to the moments, and hence the 
MC, when N is large. Recall that 5(0) = N, therefore 
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we need to consider all terms involving S (20) N p in 

(2p-l)UN p , and 

2k ~\ 



|J(0)| 2p ). We derive (|J(0)| 



-p 1 



2i 



5(20) 
N 



(4) 



Using these equations and employing the replica 
method, we obtain the log-averaged MC Lmc = 
(ln|J(0)| 2 ) - (ln|J(0)| 2 ). The typical MC of a sam- 
ple, G(0) = exp((ln | J(0)| 2 )), is then given by 



G(4>) 
G(0) 



exp(L M c) 




5(20) 



N 



(5) 



Eq. (5) is one of our main results. It provides a compact 
closed-form expression for the MC, as an explicit function 
of the magnetic flux. From Eq. (5) it becomes evident 
that a magnetic field leads to an increase in the positive 
MC: G(0)/G(O) increases from 1 to a saturated value 
2 (since 5(20) decreases from N to 0) when the flux is 
turned on and increased. G(0) = 2G(0) at the field 
that satisfies 5(20) = 0. Furthermore, it is clear that the 
MC varies periodically with the magnetic field and the 
periodicity in the flux is equal to hc/2e. 

It is illuminating to draw attention to the close 
relationship between the behaviors of lim (20) = 
5 m , m (20) and the corresponding G(0). When = 
0, (^2m(0)/iV) 2 = 1, which is the largest value of 
(/2m (20) /AO 2 as a function of 0, and the MC is equal 
to the smallest value G(0). When the magnetic field is 
increased from zero, (/2m(20)/A r ) 2 quickly approaches 
(more rapidly as m becomes larger) its smallest value, 
which is zero, at 0/2-7T = l/4m. At the same time, the 
MC rapidly increases to the largest value 2G(0). The 
physical implication of this is clear: fully constructive 
(destructive) interference in the case without disorder 
leads to the smallest (largest) hopping conduction in the 
presence of disorder. Moreover, when m (the system size 
is mxm) is large, G(0)/G(O) remains in the close vicinity 
of 2 for 0/27T > l/4m in spite of the strong very-small- 
magnitude fluctuations of /2m(20)/W around zero. 

The saturated value of the magnetic field B sa [i.e., the 
first field that makes G(0) = 2G(0)] is inversely propor- 
tional to twice the hopping length: the larger the system 
is, the smaller B sa will be. In other words, as soon as the 
system, with hopping distance r = 2m, is penetrated by 
a total flux of (l/2r) x (r/2) 2 = r/8 (in units of <J> ), the 
MC reaches the saturation value 2G(0). 

To examine the behavior of the MC in the low- 
flux limit, we first define the relative MC, AG(0) = 
[G(0) - G(0)]/G(0). Since lnJ r = -r 2 (r + l)0 2 /96 and 
ln5 r _i,i = -(r 2 - l)0 2 /24 from Eq. (3), it follows then 
that, for very small fields, in 2D AG(0) ~ \/3r 3/2 0/6 and 
in ladder-type quasi-ID structures AG(0) ~ V3r0/3. 

Our results for the MC are in good agreement with 
experimental measurements. For instance, a positive 



MC is observed in the VRH regime of both macro- 
scopic In 2 03_a; samples and compensated n-type CdSe 
||. Moreover, saturation in the MC as the field is in- 
creased is also reported in Ref. 0. 

The result for G(0) presented in this work is consis- 
tent with theoretical studies based on an independent- 
directed-path formalism and a random matrix theory of 
the transition strengths [^J. The advantages of our re- 
sults lie in that they provide (i) an explicit expression 
for the MC as a function of the magnetic field, and thus 
a straightforward determination of the period of oscilla- 
tion, (ii) explicit scaling behaviors (i.e., the dependence 
on the hopping length as well as the orientation and 
strength of the field) of the low-flux MC in quasi-ID, 2D 
and 3D systems, and (iii) allow a quantitative compari- 
son with experimental data. It is important to emphasize 
that our analytic result for the MC is equally applicable 
to any dimension, since the essential ingredient in our 
expressions is the QI quantity , which takes into ac- 
count the dimensionality. 

Quantum Interference and the small-field magnetocon- 
ductance on a 3D cubic lattice. — Let S m . n .i (= 
in 3D) be the sum over all phase factors associated 
with directed paths of m + n + I (= r) steps along 
which an electron may hop from (0,0,0) to the site 
(m,n,l) with m, n and I > 0. The vector potential of 
{B x , By, B z ) can be written as A = {zB y — yB z , xB z — 
zB x , yB x — xB y )/2. Also, a/2Tr, b/2iT and c/27r rep- 
resent the three fluxes through the respective elemen- 
tary plaquettes on the yz-, zx- and xy-planes. To com- 
pute Sm^n^i, we start from the following recursion rela- 
tion S \'l = V™ V" A i i p iO-pb)/2 c , 

where A p ^ q ^^ m ^ nt i is the sum over all directed paths 
starting from (p,q,l) and ending at (m,n,l). The phys- 
ical meaning of this relation is clear: the site (m, n, I) is 
reached by taking one step from (p,q,l — I) to (p,q,l), 
acquiring the phase i(qa — pb)/2, then traversing from 
(p, q, I) to (m, n, I) on the z — I plane. We find that 
Ap,q,l^m,n,l = exp{i[(m - p)(lb - qc) + (n - q){pc - 
la)]/2}S m -p tn - q (c). By applying the above equation I 
times, we obtain a general formula of S m ^ n ^ for m,n,l > 
1 in terms of the fluxes a, b and c as 

S m ,„,l(a, b, C) = e -i{nla+lmb+mnc)/2 jT^^ ^ c)) (6) 



C 



Pj + i Qj + i 

EE* 



,i[qja+(m-pj )b+pj (qj+l-qj)c] 



Pj + i-Pj^j + i-Qi ( C )] } ^Pl,9l( C )' CO 

with pi+i = m, qi+\ = n, and the L P! g(c)'s are de- 
fined as in Eq. (I). In the absence of the flux, S mi n t i = 
(rn + n + l)l / mini 11 = AT gives the total number of r-step 
paths connecting (0,0,0) and (m, n, /). In the very-low- 
flux limit, exactly calculated to second order in the flux 
and omitting the term In Af, the 3D analog of the har- 
monic shrinkage of the wave function becomes 
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lriiS m n i = \nla 2 + Imb 2 + nine 2 (8) 

24 

+m(lb — nc) 2 +n(mc — la) 2 + l(na — m&) 2 ] . 

In order to see how the interference patterns and the 
MC vary according to the orientation of the applied field, 
we now examine two special cases: B|| = (1, 1, l)((f>/2ir) 
and Bj^ = (1/2, 1/2, — l)((j>/2n), namely, fields paral- 
lel and perpendicular to the (1,1,1) direction. We find 
that their S m>m , m , designated respectively by xj and l^r 
(where r — 3m), exhibit quite different behaviors. Fur- 
thermore, they are insensitive to the commensurability 
of <fi, unlike the case on a square lattice. Physically, this 
can be understood because paths have a higher probabil- 
ity of crossing (and thus interfering) in 2D than in 3D; 
thus making QI effects less pronounced in 3D than in 
2D. A similar situation occurs classically (e.g., multiply- 
scattered light in a random medium). For very small 0, 
lnXjJ = -r 2 2 /72andlnX I f = -r 2 (r+l)</> 2 /144. The 3D 
behavior of AG(d>) thus become clear: ~ r</>/3 for By and 
~ \/2r 3 / 2 (f>/6 for Bj^. These results can be interpreted 
as follows: the effective area exposed to By is smaller 
(~ r), similar to our quasi-lD case with AG{<f) oc r<j>; 
while the effective area exposed to Bj^ is larger (~ r 3 / 2 ), 
thus closer to the 2D case with AG(</>) oc r 3 ^ 2 (f>. 

Average of the magnetoconductance over angles. — In 
a macroscopic sample, the conductance may be deter- 
mined by a few (i.e., more than one) critical hopping 
events. As a result of this, the observed MC of the whole 
sample should then be the average of the MC associated 
with these critical hops ||Tc(] . Thus, in 3D systems, it is 
also important to take into account the randomness of 
the angles between the hopping direction and the orien- 
tation of the applied magnetic field. We therefore con- 
sider the following picture: the ending site of all hopping 
events (with the same hopping length r) is located at 
the body diagonal (r/3, r/3, r/3) and the magnetic field 
can be adjusted between the parallel and perpendicular 
directions with respect to the vector d = (1,1,1). Our 
interest here is in the MC averaged over angles, denoted 
by AG, in the low-field limit. From Eq. (8), we obtain 
AG{B) = (27r/3%/3)rB(l + rsin 2 ^) 1 / 2 , where B is the 
magnitude of the field and u) is the angle between B and 
d. By averaging over the angle u>, we obtain 

S5(B) = ^ r VFTT BE (f,^=), (9) 

where E(ir/2, yfrfy/r + 1) is the complete elliptic inte- 
gral of the second kind. When r is large, E ~ 1 and 
we therefore have AG{B) ~ (4/3-\/3) r 3 / 2 B. This means 
that the dominant contribution to the MC stems from the 
critical hop which is perpendicular to the field. This is 
understandable through our earlier observation that the 
effective area enclosed by the electron is largest when B 
is perpendicular to d. From the above analysis, we con- 
clude that in 3D macroscopic samples the low-field MC 
should in principle behave as r 3 / 2 B. 



Finally we briefly address five issues. First, in addi- 
tion to the dominant terms, we have also obtained the 
second-order contribution to the moments (oc A^ 1 ) and 
the MC (oc l/N). The principal features in the behav- 
ior of the MC are not significantly modified: only the 
magnitude of the positive MC, including the saturation 
value, is slightly increased (though negligibly small). Sec- 
ond, we have also studied another model of disorder: 
uniformly distributed between —W/2 and W/2. The re- 
sult for the MC remains the same. Third, returns to the 
origin become important for less strongly localized elec- 
trons, and their QI effects |llj can be incorporated in our 
approach. Fourth, the main limitations of our study are 
the following: no inclusion of spin-orbit scattering effects 
(for this see, e.g., Refs. [7,8] and references therein), and 
no explicit inclusion of the correlations between crossing 
paths, as discussed in Refs. [4,7]. However, these correla- 
tions are negligible when spin-orbit scattering is present 
@. Fifth, our result for AG(B) in 3D can be tested ex- 
perimentally by measuring the MC of bulk samples for 
varying orientations of the field. 

In summary, we present an investigation of QI phe- 
nomena and the magnetic-field effects on the MC for 2D 
and 3D systems in the VRH regime. We provide exact 
and explicit closed-form results for the forward-scattering 
paths which, in the strongly localized regime, give the 
dominant contribution to the hopping conduction. 

We are very grateful to B.L. Altshuler, M. Kardar, E. 
Medina, Y. Meir, Y. Shapir, and J. Stcmbridge for their 
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